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Introduction 

The study of the fundamental theory of the strong interaction, Quantum Chromo Dy- 
namics (QCD), in extreme conditions of densities and temperatures has attracted much 
experimental and theoretical interest in recent years. Experimentally, the only avail- 
able means to investigate the many-body chromo-dynamics of a dense and hot system of 
partons involves the use of large atomic nuclei collided at relativistic energies. In this 
context, the primary goal of high-energy nucleus- nucleus collisions is the creation and 
study in the laboratory of a deconfined state of "plasma of quarks and gluons" (QGP) 
predicted by QCD calculations on the lattice [l] for values of the energy density, e > 
0.7 GeV/fm^, five times larger than those found in the nuclear ground state. 



The search of QGP signals in the A+A data relies strongly on direct comparison to 
reference results from proton-proton collisions in free space as they provide the "QCD 
vacuum" baseline to which one compares the heavy-ion results in order to extract infor- 
mation about the properties of the produced hot and dense QCD medium. Among all 
available experimental observables, hard probes provide the most direct information on 
the fundamental quark and gluon degrees of freedom. Indeed, in all hadronic collisions 
(p+P) P+A or A+A), the production of particles with high transverse momentum (jets, 
single hadrons with ^ 2 GeV/c, prompt 7) or large mass (heavy quarks), results from 
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direct parton-parton scatterings with large momentum transfer ("hard processes"). 
Since the hard cross-sections can be theoretically calculable by perturbative methods 
via the coUinear factorization theorem 0, inclusive high px hadrons, jets, direct pho- 
tons, Drell-Yan, and heavy flavors, have long been considered both experimentally and 
theoretically sensitive and well calibrated probes of the small-distance QCD phenomena. 

Simple arguments based on QCD factorization for hard cross-sections in A+A 
collisions (with its implicit premise of incoherent parton-parton scattering) [H], and 
direct experimental measurements of Drell-Yan production in Pb+Pb at CERN-SPS |3], 
and of prompt-7 jS] and total charm yields in Au+Au at BNL-RHIC, support 
that hard inclusive cross-sections in A+A reactions scale simply as times the 
corresponding hard p+p cross-sections: E da^^'^,.^x/d^P = A? ■ E do^'^_f^y^xl ^V- Foi' ^ 
given centrality bin (or impact parameter 6) in a nucleus-nucleus reaction, such a simple 
scaling rule translates into the so-called 'Wco«z (binary) scaling" relation between hard 
p+p cross-sections and A+A yields: E dNf^thxl d^P {b) = {TAA{b)) ■ ^daj^lf^^^^/d^p, 
where T^^(6) is the Glauber (eikonal) nuclear overlap function || at b. The standard 
method to quantify the effects of the medium in a given hard probe produced in a A+A 
reaction is, thus, given by the nuclear modification factor. 

"hot/dense QCD medium" d^NAA/dydpx , 

RAA{PT,y;b) = "QCD vacuum" " {TAA{b)) x d'^ajdydpr^ 

which measures the deviation of A+A at b from an incoherent superposition of A^A^ 
collisions. 



In this paper we discuss experimental results on inclusive single hadrons at high pt 
in p+p and A+A collisions at top CERN-SPS energies (\/s ~ 20 GeV) and intermediate 
RHIC energies (^i = 62.4 GeV), as well as theoretical expectations on direct photon 
production in p+p and A+A collisions at top RHIC energies (y^ = 200 GeV). The 
corresponding high px hadron and photon A+A nuclear modification factors provide, 
respectively, critical information on established QGP signals such as high pt leading 
hadron suppression (due to parton energy loss in the deconfined system) and thermal 
photon emission from the plasma. 



Case I : High pt p+p ^ 7r° + X reference at yfs ^ 20 GeV 

One of the canonical signatures of QGP formation in A+A collisions is the observation 
of a suppressed production of high pt leading hadrons (as compared to p+p collisions in 
free space) due to the non-Abelian energy loss of the parent parton traversing the dense 
medium produced in the A+A reaction ("jet quenching") [7j. In agreement with these 
expectations, high-p^ hadroproduction in central Au+Au collisions at ^/s^ = 130, 200 
GeV at RHIC has been found to be strongly suppressed (by up to a factor of 4-5) 

II The term "Ncoii scaling" comes from the fact that the average number of nucleon-nucleon (NN) 
coUisions at b is simply {Ncoii{b)) = (TAA{b)) ■ crlT^ ■ 
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[HI El Uni UH 1121 Unj compared to p+p collisions measured at the same ^/s^ [TU ITT| . 
In contrast to this result, high px pion production in central A+A at CERN-SPS en- 
ergies was found not to be suppressed but enhanced compared to production in free 
space 113 Cni EI • Such a "Cronin effect" was consistent with previous observations in 
fixed-target p+A collisions at Fermilab {y/s^ 20 - 40 GeV) [IHl CHI 1201 and in a + a 
interactions at the ISR collider {^/s^ = 31 GeV) [21] • The prevalence of the Cronin 
broadening, characteristic "cold QCD matter" systems, in A+A reactions at the SPS 
suggested that multiple initial-state (soft) kr "kicks" suffered by the colliding parton in 
their way through the nucleus dominated over possible final-state energy-loss effects for 
center-of-mass energies of order ^/s^ ~ 20 GeV. 




A recent reanalysis of the SPS results [22] has shown, however, that the apparent 
strong pion enhancement observed is not actually supported by the A+A data but 
is due to the use of inexact p+p baseline references. Indeed, since no concurrent 
measurement of high pt pion production in proton-proton collisions was carried out 
at SPS at the same ^/s as the heavy-ion experiments (and given that perturbative QCD 
calculations are not reliable in this relatively low energy domain) one had to count on 
p+p parametrizations extrapolated from data at higher collision energies. On the one 
hand, the WA98 collaboration [121 employed a empirical power-law form A \po/{pt+Po)T 
(originally proposed by Hagedorn tuned to reproduce the ISR p+p pion spectra 
(^^/s !^ 20 - 31 GeV), plus an xt scaling prescription [21] to account for the collision 
energy dependence of the cross section. On the other hand, Wang&Wang [T2| adopted 
a more complex power-law ansatz for the pt spectrum which described the charged pion 
data at y/s = 19.4 GeV combined with a pQCD parton model calculation to scale 
the cross-section down to ^/s = 17.3 GeV. Though both parametrizations were tuned to 
reproduce a subset of the existing p+p n+X data at y/s ~ 20 GeV, no true global 
analysis was carried out to fully compare the parametrizations to all the existing data 
in this energy regime. 

Figured] (left) shows all single inclusive pion spectra measured at high pt in p+p,p 
collisions in the same ^/s range of the CERN-SPS heavy-ion experiments, whereas the 
plot in Figure Q (right) shows the relevant parton kinematical domain for hard hadro- 
production in this energy regime. Fig. |2] confronts the two proposed parametrizations 
to the experimental data. Both parametrizations fail to describe adequately the shape 
of the Pt spectra, and both undershoot the cross-sections^ by as much as by factors of 
2-3 of the same order as the reported Cronin enhancements. 

As an alternative to the WA98 and Wang&Wang p+p references, we proposed [22] 
to employ a purely empirical 11-parameter functional form derived by Blattnig et. al [21] 
from a global analysis of most of the available proton-proton vr^ spectra measured in 
the range \/s ~ 7 - 63 GeV. Fig. [H] shows that the level of agreement of the Blattnig 

% Note that, as discussed in the 7r° data of [2S| plotted in Figs. n]21h£^ve been even scaled down 
by a factor of 0.8 above pT — 1.5 GeV/c to take into account the rj contamination not subtracted in 
the original tabulated results (see a discussion of this effect in Section Hi. 2. 3|) . 
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Figure 1. Left: Compilation of single inclusive pion spectra measured at high pT 
(pT ^ 2 GcV/c) in p+p collisions 1251 1^ I27j in the same collision energy range 
(yi = 16.9 - 19.4 GeV) as the CERN-SPS A+A experiments. Right: Scaling variables 
(xi^2) (average parton fractional momentum) and (z) (average momentum fraction of 
the parent parton carried by the leading pion) for p+p tt'^ production (y^ = 17.3 
GeV) at mid-rapidity versus the tt'^ momentum, computed in perturbative QCD |28| . 



parametrization to the p+p pion cross-sections of Fig. |21 is more satisfactory, both in 
shape and magnitude, than the two previous parametrizations at least within the pt 
range {pt ~ 1-5 - 3.5 GeV/c) covered by the heavy-ion experiments. Importantly, the 
p+p parametrization describes reasonable well the ^/s dependence of the yields which 
is a critical requirement for using it as a fair baseline for the Pb+Pb and Pb+Au data 
at ^/s^ = 17.3 GeV. 

Using the reference of Blattnig et. al, we revised in the nuclear modifications 
factors for the whole set of high prp data from the SPS heavy-ion programme: vr*^ and 
7r± at = 17.3 GeV from Pb+Pb (WA98) ^ and Pb+Au (CERES/NA45) ^ 

respectively, and 7r° from S+Au at ^/s^ = 19.4 GeV (WA80) [HI]. The revisited 
Raa, plotted in Figure IH indicate that high-p-r hadroproduction at ^/s^ ~ 20 GeV 
is not enhanced in central nucleus- nucleus reactions but, within errors, is consistent 
instead with scaling with the number of A^A^ collisions. Interestingly, at variance 
with central collisions, high pt pion production in peripheral Pb+Pb reactions at 
SPS is indeed enhanced by as much as by a factor of ~2 compared to the '^Ncou 
scaling" expectation ^| |221- This fact indicates that there must exist a mechanism 
in central A+A that "neutralizes" the Cronin enhancement down to values consistent 
with Raa ~ 1- Theoretical predictions jHSI of high pt tt*^ production in central Pb+Pb 
collisions including Cronin broadening and ( ant i) shadowing supplemented with final- 
state partonic energy loss in an expanding dense system with effective gluon densities 
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Figure 2. Relative differences between the single inclusive pion spectra measured 
in p+p collisions at = 16.9 - 19.4 GeV [13 EHI EHl E3 and the p+p 
TT^+X parametrizations proposed by the WA98 collaboration ^H] (upper figure) and 
Wang&Wang ^7] (lower figure) at the corresponding ^s. The shaded band represents 
the 20% overall uncertainty assigned originally to the WA98 parametrization. 
Wang&Wang only provides the fit parameters for 2 fixed values of ^/s = 17.3, 19.4 GeV. 



dN^/dy = 400 - 600 reproduce well the observed suppression factor (yellow band in 
Fig. supporting the idea that a moderate amount of jet quenching is already present 
in the most central heavy-ion reactions at SPS. The confirmation of this conclusion 
requires, however, a direct (and accurate) measurement of the high p^ p+p pion cross- 
section at ^/s = 17.3 GeV. Unfortunately, although RHIC can run Au+Au at a similar 
center-of-mass energy of ^/s^ = 19.6 GeV (corresponding to the 9 GeV injection energy 
from AGS), the minimum collision energy in the proton-proton mode is ^/s = 48.6 GeV 
(the RHIC design injection energy for each proton beam from AGS is 24 GeV, above the 
transition energy of 22 GeV) [^2] ■ This fact hinders the possibility of directly comparing 
high Pt hadropro duct ion in A+A and p+p collisions at RHIC at center-of-mass energies 
comparable to SPS. 
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Figure 3. Relative differences between the single inclusive pion spectra measured 
in p+p collisions at = 16.9 - 19.4 GeV [HI 03 EE 03 and the p+p -> tt°+X 
parametrization derived in ref. f2S|. The shaded band represents the 25% overall 
uncertainty that we assign to the parametrization. 
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Figure 4. Nuclear modification factors for pion production at CERN-SPS in central 
Pb+Pb [ini, Pb+Au IS3, and S+Au 31 reactions at ^i]^ w 20 GeV, and at ISR in 
minimum bias a + a reactions at ^/sj^ = 31 GeV The Raa from SPS are obtained 
using the p+p parametrization proposed in ref. |22| . The Raa for ISR a + a has been 
obtained using the p+p spectrum measured in the same experiment [IJ. The shaded 
band at Raa — 1 represents the overall fractional uncertainty (including in quadrature 
the 25% uncertainty of the p+p reference and the 10% error of the Glauber calculation 
of Ncou)- CERES data [21] have an additional overall uncertainty of ±15% not shown 
in the plot. 
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Figure 5. Nuclear modification factors for pion production at CERN-SPS in central 
Pb+Pb [ni, Pb+Au inni, and S+Au [JT] reactions (Fig.|3I) compared to a theoretical 
prediction from Vitev and Gyulassy |321 including standard nuclear effects (Cronin 
broadening and shadowing) and final-state parton energy loss in a system with initial 
gluon densities dN^/dy = 400 - 600. 



Case II : High px p+p vr°, h"^ + X references at y/s = 62.4 GeV. 

The study of the excitation function of high px hadron suppression between top SPS 
and top RHIC energies was the main motivation behind the dedicated Au+Au run at 
RHIC intermediate energies {^/s^ = 62.4 GeV) carried out in April 2004. PHENIX 
measured charged hadrons up to pr = 5 GeV/c and identified neutral pions in the 
range pt = 1 - 7 GeV/c [211 • PHOBOS [23 and STAR jSni measured inclusive charged 
hadrons up to pt ~ 4.5 GeV/c and 12 GeV/c respectively. However, as in the SPS case, 
no concurrent p+p reference measurement was performed at ^/s = 62.4 GeV, and the 
corresponding Au+Au nuclear modifications factors were constructed using, basically, 
p+p h^,TT^ differential cross-sections measured in the 70s and 80s at the top ISR 
energies {^/s = 62 — 63 GeV). Tabled collects all the existing measurements of highp-r 
(neutral and charged) pion and inclusive charged hadrons in proton-proton collisions at 
the maximum CERN-ISR energy. 

II. 1. p+p — s> /i^ + X reference at ^/s = 62.4 GeV 

Rather than using the existing ISR data at ^/s = 62 — 63 GeV, the PHENIX p+p in- 
clusive charged hadron reference at a/s = 62.4 GeV has been independently obtained [S2] 
by interpolating from lower and higher collision energy measurements'*': p+p at ^/s = 

+ The same procedure was followed to obtain the p+p reference for RHIC Run-1 Au+Au measurement 
at = 130 GeV 0. 
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Table 1. Chronological compilation of 7r°'* mP = tt" + ?]], and inclusive charged {h^) production measured at the top CERN-ISR 
energies. For each (1) reaction, we quote the (2) center-of-mass energy, (3) collaboration/experiment name, (4) bibliographical reference, 
(5) rapidity domain, (6) measured pt range (center of min. and max. bins quoted), (7) total number of data points, (8) tt'^ energy-scale 
errors (the value in parenthesis is the true energy-scale uncertainty and the uncertainty quoted is the total effect of the error propagated 
to the yields), (9) additional systematic and/or absolute normalization (luminosity) errors [(*) Kourkoumelis et al. energy scale error 
includes all syst. uncertainties and is an average over pr (error values quoted in the paper are in the range 27%-42%). (**) Akesson et al. 
original point-to-point errors include all uncertainties.]. Column (10) lists those measurements that have been corrected to account for the 
(slightly) different ISR {y/s = 62. - 63. GeV) and RHIC {^/s = 62.4 GeV) center-of-mass energies (those data sets marked with a dagger 
{^) have not been revised since they cover low pr values where the effect of the correction on the yields is minimal). Columns (11) and 
(12) indicate whether the published 7r° spectrum has been corrected for direct-7 and ij "contaminations" resp. as described in the text. 
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21, 31, 44, and 53 GeV at CERN-ISR gZI, PHENIX p+p at 200 GeV, and CERN-UAl 
p+p data at 200, 500, and 900 GeV [HZ|. Fig. El shows an example of the cross-section 
interpolation procedure in four individual pt bins from the lower and higher ^/s mea- 
surements. The resulting interpolated spectrum at 62.4 GeV is fitted to a modified 
power-law form: Ed^app-^hx/d^P = A{e"''^ +Pt/po)~"; with ]9re/immar?/ parameters A 
= 196.4 [mb GeV-^c^], a = 0.0226 [GeV-^c], po = 2.301 [GeV/c], and n = 14.86; with 
an assigned overall uncertainty of ±25% [22] ■ As an independent cross-check the fit is 
compared (Fig. to the three ISR /i^ data sets measured at 63 GeV (three bottom 
rows of Table Q). The agreement data-fit is good within the assigned ±25% uncertainty 
of the parametrization. 




Figure 6. Charged hadron cross-sections as function of center of mass energy for fixed 
Pt values from 2.5 to 4 GeV/c |52]- Tlie lower energy points (red) are from fits to ISR 
data (combined tt^, K and proton spectra from Alper et al. |47j). the higher energy 
points (blue) are from fits to UAl data [S7|, and the black points are from a fit to 
PHENIX data. 

PHOBOS fitted the experimental p+p data from the CDHW experiment jSI] 
at ?7 = 0.75 (the same rapidity range of their spectrometer) to the following expression: 
Ed^app^hx/d^p = A{1 + Pt / Po^ Pt / \/ pI + a, with A = 244.5 [mb GeV-^c^], po = 
2.188 [GeV/c], a = 0.0085 [GeV^/c^], and n = 15.37. STAR preliminary lefeieiace p+p 
spectrum jSEj used a Hagedorn power-law form, A[l +Pt/po]~"', with parameters A 
= 292.48 [mb GeV-^c^], po = 1.75 [GeV/c], and n = 13.23. Fig. [7| compares the 3 
parametrizations to the existing ISR inclusive charged hadron cross-sections at 63 GeV. 
The 3 parametrizations agree among each other (and the data) within ±20%. A direct 
measurement of the inclusive charged hadron production at large pt in p+p collisions 
in a dedicated run at y/s = 62.4 GeV at RHIC is mandatory, however, if one wants to 
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reduce the corresponding uncertainties propagated to the Au+Au nuclear modification 
factors and constrain more quantitatively the model predictions for the parton energy 
loss excitation function (see Section ITL3|1 . 
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Figure 7. Upper: Inclusive charged hadron spectra measured in p+p collisions at -y/s 
= 63 GeV at CERN-ISR [SUIEII at midrapidity {g = 0, 0.25) and moderately forward 
rapidities (77 = 0.75) compared to PHENIX |S2|, PHOBOS and STAR ^ /i± 
parametrizations. Lower: Ratios of the same experimental spectra over the three 
parametrizations (each one evaluated in their valid range of rapidities). 
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II. 2. p+p 



+ X reference at -^i = 62.4 GeV 



In order to obtain a benchmark p+p— tt" + X reference spectrum at ^/s = 62.4 GeV we 
first collected (Table Hj) all existing experimental 7r° (9 measurements pTll^l^HUll^ 
iailSlllllliaiini), tt^ (4 measurements gOlllTlllHllini) and (3 measurements EDIEII) 
at the highest ISR collider energies {^/s = 62 - 63 GeV) and added in quadrature the 
original systematic and normalization uncertainties to the point-to-point errors. We 
included in our compilation the averaged (7r+ + n~)/2 spectra* as well as the inclusive 
charged hadron spectra divided by the measured /i^/vr = 1.6 ± 0.16 ratio (Fig. |HI) since 
this provided us with additional constraints for our global fit analysis at relatively low 
Pt where no neutral pion data is available. The corresponding data points (adding to 
a total of ~300) were fitted to a common functional form. The ratio of the weighted 
average fit over each data set is shown in the upper plot of Fig. IHl Large differences in 
the shape of the px distributions and in the magnitude of the cross-sections are evident 
which implies that many of the measurements are inconsistent among each other well 
beyond the originally quoted uncertainties. 




4 5 
(GeV/c) 



Figure 8. Ratio of total charged hadron over pion spectra in p+p colhsions at ISR 
energies (y/s = 21, 31, 44, 53, and 63 GeV The straight line is at h^/w = 1.6. 

Investigation of the original published results indicates, however, several relevant 
differences affecting the experimental measurements. First, although at first sight 
discrepancies of order ~1 GeV in the center-of-mass energy should not dramatically 
modify the single 7r° spectra, it turns out that at very high pj- the absolute difference 
in the perturbative yields between ^/s = 62 GeV and 63 GeV can indeed be as large as 
~20% fFig.lTUj). Secondly and most important, only one experiment (Akesson et al. [IHj) 
fully identified neutral pions via a standard invariant mass analysis of photon pairs. The 
rest of the experiments either did not separate the 77 decay of 7r° and rj and/or did not 

* Isospin symmetry justifies the assumption: tt*^ k. (tt^ + 7r^)/2. 
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subtract the direct photon component from the experimentally measured "unresolved" 
vr" spectra. In order to do a meaningful comparison and average of all 7r° data sets, one 
needed therefore to subtract each one of these experimental "contaminations" from the 
published data tables. 
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Figure 9. Ratio of all experimental measurements reported in Tabled over the final 
p+p 7r° parametrization (Eq. |2Jl). The upper plot shows the data as originally 
published whereas the lower panel shows the corrected data as described in the text. 
The error bars for each data set include (in quadrature) the original point-to-point, 
systematic and normalization uncertainties. 
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II. 2.1. Center- of -mass energy correction: CERN ISR was a collider with proton beams 
crossing (and colliding) with a non-null angle depending on the experimental runs and 
setups. This fact resulted in differences in the effective available energy in the center- 
of-mass as large as ~1 GeV (from ^/s = 62 to 63 GeV). RHIC Au+Au collisions were 
instead performed at a fixed y^s^ = 62.4 GeV. Fig. [101 shows the perturbative ratio 
of 7r° cross-sections as a function of pr for p+p collisions at = 62 and 63 GeV over 
those at ^/s = 62.4 GeV as given by NLO calculations Whereas the effect on 

the yields is minimal (^5%) for p^ < 8 GeV/c, the difference in yields monotonically 
increases with pt reaching a maximum of ~10% at the highest measured pt values at 
ISR. We corrected the yields of all data sets indicated in the last column of Table [T] 
using a simple second-order polynomial pt fit of the computed pQCD ratio. 




Figure 10. Perturbative ratios of p+p 7r° yields vs. pr at two different center-of- 
mass energies (y/s = 63 GeV, upper curve, and 62 GeV, lower curve) over the yields at 
y/s = 62.4 GeV, as given by NLO pQCD calculations (note that any uncertainties 
in the PDFs and/or FFs basically cancel out in the ratio). 



II. 2. 2. Direct photon subtraction: Direct 7 were discovered in p+p at CERN-ISR 
in 1979 54,:, therefore before this date any high-p-p photon-like cluster detected in 
the electromagnetic calorimeters was "identified" as a neutral pion. Figure ^2 shows 
the direct-7/7r° ratio as a function of pt measured at y/s ~ 62.4 GeV by three ISR 
experiments (squares) [SH ESI US] compared to the NLO pQCD ratio (circles) computed 
with CTEQ6 PDF and for three different (factorization-renormalization) scales (the 
theoretical points are centered at fi = pr and the "errors" correspond to fi = pt/ 2 - 
2pt) I56j . The prompt photon "contamination" is marginal below ~4 GeV/c, but 
it accounts for ~l/3 of the vr" yield at pt ~10 GeV/c, and equals it at ~14 GeV/c. 
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The agreement data-theory is good and allows us to extrapolate the ratio to px values 
higher than those measured in |3U ESJ 03] . The combined experimental and theoretical 
data points have been thus fitted to a 4th order polynomial (black curve) as a function 
of pt, R^/Troipr), with parameters: po = 4.55e-02, pi = -6.04e-02, p2 = 2.51e-02, p3 = 
-2.53e-03 and ps = 1.03e-04. The corrected 7r° yields, Yt,o{pt), were obtained from the 
"unresolved" 7r° yields, Yj,o^^{pt), via: F^o = ¥^^0+^ ■ R~J^o/(l + R~}_n), for those data 



sets listed in the 11th column of Table [T] 
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Figure 11. Ratio of j/tt'^ cross-sections in p+p collisions at y/s « 62.4 GeV as a 
function of pt measured experimentally |54l I55L 1^ (squares), and computed in NLO 
pQCD (circles, with error bars covering the range of theoretical scale uncertainties: ji 
= Pt/2 — "^Pt) inSlinni' The black curve is a fit of the experimental and theoretical 
results to a common 4th order polynomial. 



II. 2. 3. — > 77 subtraction: Many of the "unresolved" 7r° measurements at ISR 
assumed that all detected electromagnetic clusters at high pt were (merged) photons 
from the 7r° decay and neglected any possible contribution from the 2-gamma decay 
channel of the rj meson (BR^^^.^ = 0.394). At high pt {pt > 1-5 GeV/c), however, the 
77/71° ratio is Rrj/nO = 0.46 (see the "world" compilation in Fig. ll2|) . Those "unresolved" 
TT^ spectra listed in (the 12th column of) Table ^ have been scaled down by a factor 
of 0.82 above pt = 1.5 GeV/c to take into account the (BR^^^^ ■ -R^/vrO = 0.394 ■ 0.45 
= 0.18) rj contamination [Note that such ~18% factor was mentioned in some of the 
original papers (e.g. jHZllSH]) but not actually subtracted from the tabulated results]. 
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Figure 12. Collected systematics of 7//7r° ratios in hadronic collisions. The average 
above pr — 1-5 GeV/c is i?^/7r0 — 0.46. 



II. 2. 4- Final parametrization: After correction of the experimental spectra as described 
in the previous sections, there remained still a few outliers measurements that were 
(partially or totally) inconsistent (beyond ~ 1.5cr) with the rest of the data (see lower 
panel of Fig. IHl). Three data sets: Drijard et al. jlH] (charged pions measured off central 
rapidity, at 50") , Eggert et al. jSHl and Busser et al. |SZI were excluded of the final global 
fit analysis. Fig. shows all the (corrected) pion data sets plotted as a function of pt 
compared to a common (purely empirical) 5-parameter functional form: 

E£app^^x/d'p = A (e'^-PT+b-PT + p^/p^)-, (2) 

with parameters: A = 265.1 [mb GeV-^c^], a = -0.0129 [GeV-^c], b = 0.04975 [GeV/c], 
Po = 2.639 [GeV/c], and n = 17.95. Eq. Q provides a very good reproductionjl of the 
full spectral shape in the range pt = - 16 GeV/c. Above pt ^ 8 GeV/c the data 
(and the fit) have an exponential-like shape. This departure from the pure power-law 
behaviour expected for parton-parton scatterings is due to the fact that such large pt 
values are in a kinematical domain above (xt) ~ 0.3 and (z) ^ 0.8 (see Fig. ^[ right) 
where both the parton distribution functions (PDFs) and fragmentation functions (FFs) 
resp. start to decrease due to phase space boundaries (PDFs and FFs reach zero at the 
kinematical limit: pt ~ 30 GeV/c). Such a change in the pt shape is confirmed also by 
pQCD calculations (see Fig. IT^. 

Figure ITU shows the ratio of the (selected and corrected) experimental pion spectra 
over the final p-|-p — 7r° + X parametrization in the same pt range {pt = 1-8 GeV/c) 

jj Mind that this is a purely empirical fit whose validity beyond the considered pt range is not 
guaranteed. 
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Figure 13. Compilation of all pion transverse spectra measured in p+p collisions at 
-^s sa 62.4 GeV and fitted to a common function, Eq. (0), with the parameters quoted 
in the text. 

covered by the PHENIX Au+Au — > 7r° + X data at RHIC. All p+p cross-sections are 
consistent with the final fit within their associated errors. The dashed lines indicate 
the ib25% systematic uncertainty assigned to the reference. The gray thick line is the 
PHENIX charged-hadron reference (see Section FlLlj) divided by the expected h^/n^ = 
1.6 ± 0.16 ratio at px > 1.5 GeV/c (see Fig.|Hl)- The ±0.16 errors of the h^/ir ratio are 
shown as dashed gray lines. The additional systematic 25% error in the reference is 
not plotted. 
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Figure 14. Ratio of the selected and corrected tt"'* and (scaled) /i measurements at 
CERN-ISR over the final 7r° parametrization (Eq. ^ with fit parameters reported in 
the text) as a function of pt- The gray line is the charged-hadron reference spectrum 
(Section Uni, divided by h^/i: = 1.6 ±0.16 (The errors of the /i^/tt ratio are shown 
as dashed gray lines. There is an additional systematic 25% error in the /i^ reference 
not plotted). 



Figures fT5llT6l show a comparison of the empirical vr'' parametrization to NLO 
pQCD predictions from W. Vogelsang [HS] with fixed PDFs (CTEQ6), 2 different sets 
of fragmentation functions (Kniehl-Kramer-Potter KKP [SH] and Kretzer |59J and 3 
(factorization-renormalization-fragmentation) scales. The general shape and overall 
magnitude of the px spectrum is well reproduced by the theoretical calculations (the fit 
is well contained within the theoretical limits given by /i = pt/2 — 2]?^)- KKP FFs seem 
to reproduce better the magnitude of the cross-section, and the scale n = pt provides 
the best agreement with the data, especially above ~ 5 GeV/c. Below ^ 5 
GeV/c all NLO spectra tend to consistently underpredict the observed 7r° cross-section. 
Such a discrepancy is also observed in high px hadro-production at lower ^/s [^01 where 
soft-gluon resummation corrections jHI] and additional non-perturbative effects (e.g. 
intrinsic kx must be introduced to bring parton model analyses into accord with 
data. 
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Figure 15. Left: Comparison of the final empirical p+p ^ tt'^ + X parametrization 
at = 62.4 GeV (solid curve) to NLO pQCD calculations for 2 sets of FFs 
(KKP 58 and Kretzer 59 ) and 3 different scales fi = pt/2,pt,2pt- Right: Scaling 
variables {xi^2) (average parton fractional momentum) and (z) (average momentum 
fraction of the parent parton carried by the leading pion) for p+p — > tt*^ (^/s = 62.4 
GeV) at mid-rapidity versus the 7r° momentum, computed in perturbative QCD [2B1- 
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Figure 16. Ratio of the NLO pQCD predictions for inclusive p+p n'^ + X 
production at ^/s — 62.4 GeV (Fig. I15|) over the final empirical 7r° parametrization 
(Eq. ||21with the fit parameters reported in the text) as a function of pr- 
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II. 3. Nuclear modification factor at ^/s^ = 62.4 GeV: 

Figure IT7I shows the preliminary nuclear modification factor, Eq. for high px 7r° 
production in 0-10% central Au+Au collisions at ^/s^ = 62.4 GeV jSlI obtained 
using: (i) the p+p vr° reference Eq. (red circles), and (ii) the p+p h^ reference 
(see Section lILlll divided by the expected h^/n^ = 1.6 ratio (open circles). [For 
complementary info, the average Glauber number of A^A^ collisions is (Ncoii) = 845.4 
±140 for a p+p inelastic cross-section at a/s = 62.4 GeV of ainei = 35.6 ± 0.5 mb 
obtained from the weighted average values of the total jHSl HH IMl tHHl IHHl EZl inil (ctot = 
43.37 ±0.17 mb) and elastic [03113 (c"ez = 7.75 ±0.10 mb) cross-sections measured at 
CERN-ISR]. Both Raa are compared to parton energy loss predictions for the quenching 
factor in a system with effective gluon densities dN^ /dy = 650 - 800 (yellow band) |70j . 
At intermediate ~ 2 - 5 GeV the use of one or the other parametrization results in 
differences as large as ~25% in the amount of the suppression. Those divergences are 
indicative of the systematic uncertainty of the obtained p+p baseline references at ^/s 
= 62.4 GeV. Clearly a dedicated RHIC proton-proton run at this collision energy would 
help to reduce these uncertainties and better constraint the theoretical predictions of 
the excitation function of the high pt suppression lU 113 IZS| • 
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• PHENIX Preliminary (0-10%) [p+p tiO ref.] 
O PHENIX Preliminary (0-10%) [p-i-p h*/1 .6 ref.] 
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Figure 17. Preliminary PHENIX nuclear modification factor, Raa{pt), ioi tt" 
measured in central Au-|-Au at 62.4 GeV Pl] obtained using the p+p tt'^ + X (red 
circles) and p+p + X (open circles) references discussed in the text, compared 

to theoretical predictions for parton energy loss in a dense medium with dN^ /dy = 
650 - 800 JO,. 
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Case III : p+p — 7 + X reference at y/s = 200 GeV 

Thermal (real or virtual) photons emitted in high-energy A+A reactions provide direct 
information on the thermodynamical properties of the radiating underlying QCD matter 
and have long been considered privileged signatures of QGP formation [21]. Direct 
photons, defined as real photons not originating from the decay of final hadrons, are 
emitted at various stages of a A+A reaction. Three qualitatively different mechanisms 
are usually considered: (i) prompt 7 ( "pre-equilibrium" or "pQCD") emission from 
perturbative parton-parton scatterings in the first tenths of fm/c of the reaction, and 
(ii) subsequent emission from the thermalized partonic (QGP) and (iii) hadronic (hadron 
gas, HQ) phases of the reaction. The partonic diagrams contributing in leading-order 
to photon production are g(?-Compton and gg-annihilation (Fig. ITHj) and coUinear g, g 
fragmentation (Fig. ITTHl. 




Figure 18. Compton and annihilation diagrams for direct photon production in 
parton-parton scatterings. 




3 <QSiSSlSSSiS>&^ 



Figure 19. Fragmentation (or "bremsstrahlung" ) diagrams for direct photon 
production in parton-parton scatterings. 

The expected photon spectrum from Au+Au reactions at ^Js^n = 200 GeV, can 
thus be obtained theoretically by combining: (i) NLO pQCD calculations for the 
primordial (hard) production (perturbative p+p yields [Sn] scaled by the nuclear overlap 
function Taa), plus (ii) hydrodynamical calculations of the space-time evolution of the 
reaction [751 US EZ] complemented with modern parametrizations of the QGP [TSj and 
HG [Zni photon emission rates. Such calculations indicate that QGP thermal emission 
should be visible in a window of the inclusive direct photon spectrum between p^ ^ 1 
- 3 GeV/c (Fig.EOI). 
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Figure 20. Expected thermal and prompt photon spectrum for central Au+Au 
reactions at ^/snn = 200 GeV as given by a hydrodynamical model calculation |77] 
complemented with pQCD yields for the prompt 7 '56'. 



However, unfortunately, this range of transverse momenta potentially presents 
difficulties for the extraction of the thermal component. Indeed, perturbative QCD 
calculations indicate (Fig. l2T| that below pt ^ 3 GeV/c, prompt photons are mainly 
produced via the parton bremsstrahlung mechanism ("anomalous" component) [At 
increasingly highp-r, prompt photon production is dominated by the purely perturbative 
production mechanisms (Compton and annihilation. Fig. I18j) but still ~l/4 of the total 
7's appear to come from jet fragmentation (diagrams depicted in Fig. IT^ in these 
calculations]. Such a component should be depleted in central Au+Au at ^/s^ = 
200 GeV due to the same final-state QCD medium effects that result in the observed 
factor of ~4-5 suppression of high pt hadroproduction. To get a handle on the possible 
effect of parton energy loss in the total prompt photon spectrum we plot in Fig. 1221 the 
"photon nuclear modification factor", R]^^{pt), for central Au+Au estimated simply 
assuming that the suppression factor for the 7-fragmentation component is the same as 
that observed for high px hadrons, i.e. R'^JX'" ~ Raa'"^ ~ 0.25. We determine i?^^ 
= [Eda-y^^Jdp — Q.7b Eda^j^^J dp\/[Eda^^^J dp] with the same NLO yields jSE] used in 
Figure EU 

Figure 1221 indicates that one should expect a moderate ~30% effective suppression 
of the total inclusive high pt photon spectra due to parton energy losses in the dense 
mediumff. This estimate is consistent with more involved calculations of the same 

ffNote that in this estimate we have not considered partially counteracting (small) effects such as 
"Cronin enhancement" and Au PDF shadowing effects [HOI- 



Relevance of baseline hard p+p spectra for high- energy A+A physics 



22 




Figure 21. Relative contribution of different subprocesses of direct 7 production 
versus pt in p+p collisions at y/s — 200 GeV according to NLO pQCD pHj : 
Compton+Annihilation diagrams (upper band) and fragmentation diagrams (lower 
band). 
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Figure 22. Nuclear modification factor, i?^^(pT), for direct 7 in central Au+Au 
collisions at ^/Jnn = 200 GeV assuming the same quenching factor for the jet- 
fragmentation photon component as observed for high px hadrons {Raa ~ 0.25). 
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effect jHIlE21- The suppression of the jet bremsstrahlung component could, therefore, 
partially "mask" the enhancement due to thermal photon emission in the range pt 
= 1-3 GeV/c. Since the jet-fragmentation 7 component cannot be experimentally 
discarded via the standard "isolation" method due to the large soft background in 
Au+Au collisions, the only way to disentangle experimentally the counterbalancing 
effects of the thermal and quenched prompt 7 in Au+Au requires a detailed analysis of 
the p+p reference: 

(i) Measurement of isolated photons in p+p, Npp^^.^^^^^^^, down to ~ 1 GeV/c 
with uncertainties ^15% provides a handle on the actual Compton+Annihilation 
reference production. 

(ii) Measurement of total inclusive photons in p+p, Npp^^^^^^^, down to ~ 1 GeV/c 
with uncertainties ^15% provides a handle (since Npp^^^^^ = Npp^^^^^^+ Npp^^^^^^^J 
on the actual fragmentation 7 reference production. 

(iii) Measurement of total inclusive 7 production in Au+Au, NAuAu^'ytotan down to px ~ 
1 GeV/c with uncertainties <15%. 

(iv) The upper limit on thermal production for a given Au+Au centrality (with nuclear 
overlap function Taa) is given by: NZau^^,^,,^^, = NauAu-.^,,,^,-Taa- Npp_.^^^^^^^^^. 

(v) The lower limit on thermal production for a given Au+Au centrality (with nuclear 
overlap function Taa) is given by: NZau-.^,^,,.^,^ = NauAu^^,,,,^ - Taa ■ Npp^^^^^^^. 

Figure 1221 shows preliminary measurements by PHENIX of the tota/ prompt photon 
production in p+p 83^ [^pp^itotai item (ii)] and in Au+Au [H] [N AuAu^-ytotai item 
(iii)] collisions at RHIC obtained by statistically subtracting the hadron decay-photon 
contributions (vr*^, rj, ...) from the total measured 7 spectrum. Within uncertainties, 
both results are consistent with the perturbative QCD expectations [56 . A separation 
of the p+p fragmentation-7 component is under- way too jHS]- The main issue in order 
to experimentally resolve a possible thermal component in Au+Au is to measure with 
small uncertainties the prompt photon production below ~ 3 GeV/c (where the 
background of decay photons and (anti)baryon contaminations is more significant) in 
both colliding systems. Complementary methods to the statistical subtraction one (e.g. 
the measurement via 7 "conversion"), that are efficient for direct photon identification 
in the interesting range px ~ 0.5 - 2 GeV/c, are being currently pursued. 

Summary 

Three different cases of hard production in proton-proton collisions at three different 
center-of-mass energies, ^/s ~ 20, 62 and 200 GeV, have been discussed as benchmarks 
for the investigation of QGP signals in nucleus-nucleus collisions. High px neutral pion 
production in central A+A collisions at CERN-SPS ( ^Smn ~ 20 GeV) has been found 
to be (slightly) suppressed compared to a parametrization of the vr*^ cross-sections in 
p+p collisions in free space. Such a result, possibly indicative of jet quenching effects 
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Figure 23. Left: Preliminary p+p — 7 + X measured by PHENIX at ^/s = 200 
GeV [S3]. Right: Preliminary direct-7 excess measured by PHENIX at ^/s]^ = 200 
GeV |2| ■ Both results are compared to NLO pQCD calculations ^56, . 

at these energies, should be confirmed by a direct experimental measurement of high pt 
p+p hadropro duct ion at ^/s ^ 20 GeV. Secondly, high pt vr° and h^ reference spectra 
at i/i = 62.4 GeV have been constructed based on the weighted averaged of revised 
experimental data from CERN-ISR. Hard hadro-production is seen to be suppressed by 
up to a factor of ~3 in central Au+Au collisions at RHIC at these collision energies. 
However, the current p+p references have px dependent uncertainties of order ~25% 
precluding a detailed quantitative study of the excitation function of high pt nucleus- 
nucleus suppression between SPS and RHIC energies. A dedicated RHIC p+p run at 
y/s = 62 A GeV would be needed to reduce such uncertainties and better constraint 
the theoretical models of parton energy loss in dense QCD matter. Finally, we have 
studied the possible effects of parton energy loss in the direct photon contribution of 
jet-fragmentation origin in Au+Au collisions at ^/s^ = 200 GeV. Estimates based on 
NLO pQCD indicate that the nuclear modification factor can be reduced by ~30% due 
to these effects, partially hiding the expected thermal photon emission from a radiating 
QGP. We have discussed how to get a handle on a possible thermal radiation Au+Au 
signal by detailed measurements of the isolated and non-isolated direct photon spectra 
in baseline p+p collisions at = 200 GeV. 
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